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ABSTRACT 

We use observations of the C I, C II, H I, and H2 column densities along lines 
of sight in the Galactic plane to determine the formation rate of H2 on grains and 
to determine chemical reaction rates with Polycyclic Aromatic Hydrocarbons. 
Photodissociation region models are used to find the best fit parameters to the 
observed columns. We find the H2 formation rate on grains has a low rate (R ~ 
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1 x 10 -17 cm 3 s _1 ) along lines of sight with low column density (Ay ~ 0.25) and 
low molecular fraction (/h 2 ~ 10~ 4 ). At higher column densities (0.25 < Ay < 
2.13), we find a rate of R ~ 3.5 x 10~ 17 cm 3 s _1 . The lower rate at low column 
densities could be the result of grain processing by interstellar shocks which may 
deplete the grain surface area or process the sites of H + H formation, thereby 
inhibiting H2 production. Alternatively, the formation rate may be normal, and 
the low molecular fraction may be the result of lines of sight which graze larger 
clouds. Such lines of sight would have a reduced H 2 self-shielding compared to 
the line-of-sight column. We find the reaction C + + PAH~ — > C + PAH is best 
fit with a rate 2.4 x IO^fahT^ 05 cm 3 s" 1 with T 2 = T/100 K and the reaction 
C + + PAH -> C + PAH+ is best fit with a rate 8.8 x 10^ 9 $pah cm 3 s -1 . In high 
column density gas we find $pah ~ 0.4. In low column density gas, $pah is less 
well constrained with $pah ~ 0.2 — 0.4. 

Subject headings: astrochemistry-ISM: clouds-ISM: general — ISM: molecules 



INTRODUCTION 



Much of the chemistry in the interstellar medium (ISM) proceeds via two body gas phase 
reactions; however, based on reaction time arguments it is well established that the formation 
of H 2 must proceed predominantly by grain surface reactions (IHollenbach &: Salpeterlll97ll ) 
and thus the equilibrium abundance of H 2 in diffuse gas depends on the balan ce between 



dissoc iation from ultraviolet radiation and the formation on grains. In addition, iLepp et al. 



( 119881 ) suggested that in diffuse clouds, reactions with small grains or large molecules can 



affect the charge balance and abundance o 



Bakes &: Tielens 



2003 



1998 



Weingartner fc Draine 



metal ions. More recent theoretical 



2001 



Chang et al.l 120061 ) and laboratory, (ICazaux fc Tielensl 120041 ) investigations have ex- 



Wolfire et al.l 20031 : iLipshtat fc Biham 



plored the rates of H 2 formation and chemical reactions on surfaces of small grains and 
their dependence on environmental factors. The chemical rates have been difficult to pin 
down in part because the exact abundance, type, surface area, and physical properties of the 
grains/large molecules are not well determined. 

Studies with a more observa tional approach have also estimat e d H 2 formation ra te s and 
chemi cal rates with grains (e.g., I Jural Il975l : iBrowning et al.ll2003l ; IWelty et al.l 120031 ) . lJura 
( 119751 ) estimated a formation rate of about R ~ 3 x 10 -17 c m 3 s~ l in diffu s e gas of low column 
density; however, one could see in the Copernicus survey of ISavage et al.l (119771 ) that a single 
H 2 formation rate can not explain all of the observed column densities even accounting for 
the expected variation in ultraviolet radiation field and gas density. A more recent analysis 
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by iGry et al.l (120021 ) using data from the Far Ultraviolet Spectros copic Explorer (FUSE) 
find a formation rate of R ~ 4 x 10~ 17 cm 3 s _1 in diffuse gas and Browning et al.l (120031 ) 
estimates the rate at 2.5 - 1 times lower i n the l ow metallicity environments of the Large and 
Small Magellanic Clouds. lHabart et al.l (120041 ) suggests that the H 2 formation rate varies 
in molecular clouds exposed to intense radiation by a factor of about fi ve, perhaps due to 
the t emperature dependence of the formation process on grain surfaces ( iCazaux fc Tielens 
20041 ) . In this paper we take a semi-empirical approach by allowing the ensemble of recent 
observations to guide us in the proper "astrophysical" rates rather than use rates based on 
laboratory measurements with an assumed grain type and distribution. 



Bakes fc Tielensl ( 119981 ) investigated the effects of Polycyclic Aromatic Hydrocarbons 
(PAHs) on the chemistry of dense regions exposed to intense ultraviolet radiation. They 
found that charge exchange reactions with PAHs are imp ortant in setting the column den- 
sity of neutral metals, in particular C I. IWelty et al.l (120031 ) found patterns of neutral column 
densities in diffuse clouds higher than suggested by electron recombination alone and sug- 
gested that recombination on grains might be important. IWeingartner fc Draind (120011 ) 
carried out a detailed analysis of the "grain assisted" rates require d to explain the observed 
column densities along the line of sight to 23 Ori. IWolfire et al.l ( 119951 . 120031 ) found that 
reactions with PAHs can affect the ionization balance of metals and also atomic hydrogen, 
in diffuse gas, in particular in the warm neutral medium. This has an important effect on the 
grain charge and thus the heating produced by the grain photoelectric effect; the dominant 
heating process in the diffuse ISM. 



In lWolfire et al.l (120031 ) we adopted rates with PAHs based on the lDraine fc Sutinl (119871 ) 
formalism for interactions with small grains, but modified the rates to match a preliminary 
assessment of the observed N C i/N cn ratio in diffuse gas. This ratio is a particularly good 
diagnostic because carbon is abundant, the rates of C + recombination on PAH anions effec- 
tively compete with electron recombination, and carbon exists in the diffuse ISM outside of 
H II regions as either C I or C II but not higher ionization stages. 



For several ye ars, t he Copernicus column densities of H I (iBohlin et a. 



1978), H 2 



(jSavage et al.lll977f) . C I (I Jenkins &: Shayalll979l ; I Jenkins et al.lll983l ). and C II ( IHobbs et al. 
19821 ) were the only data available. More recently the combined observations from the Hub- 
ble Space Telescope Goddard High Resolution Spectrograph (GHRS) and Space Telescope 
Imaging Spectrograph (STIS) instruments along with the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE) have produced many mor e sight li nes to study in particular in H? (e.g., 



Rachford et al. 2002 



Sofia et al. 



Cartledge et al.ll2004ri . C I (e.g.. lZsargo et al.lll997l : IZsargo fc Federman 



2003h . and C II (e.g 

and in § 3 we discuss our cloud models and our fitting procedure 



200J). In 



2 we present the observations used in this study 

In § 4 we present our 
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results for the best fit rates for grain reactions and for H2 formation. We also discuss some 
of the interesting sight lines which seem to have outlying points. We summarize our main 
conclusions in 5 5. 



OBSERVATIONS 



Observations of the H I, H 2 , C I, and C II columns are listed in Tableland are taken 
from the literature as noted. We have included only those lines of sight which have measured 
values of logiVni, log A^h 2 , and logiVci, and have not included lines of sight with upper or 
lower li mits in these quan tities. The C II column is measured directly only towards thirteen 
sources ( jSofia et al.l 120041 ) . Six of these have definite C I columns and are included in our 
sample. For all other sources we use the conversion Nqii = Ac x N where N = Nm + 2N# 2 , 
and .Ac = 1.6 x 10~ 4 is the mean gas phase abundance of carbon per hydrogen nucleus found 



by ISofia et al.l (120041 ) from all measured values of the C II carbon abundance in the diffuse 
ISM. 

Most of the H I and H2 columns a nd uncertainties are from the Copernicus ult raviolet 
absorpti on spectroscopy o bserv ations of iBohlin et al.l (119781 ) and lSavage et al.l (119771 ). More 
recently, ICartledge et al.l (120041 ) updated several of these lines of sight with FUSE observa- 
tions. The columns are in good agreement wi th the previous results, but with smaller quoted 
uncertainties. In cases of overlap we use the ICartledge et al . (2004) results. Additional H 2 
columns are provide d by FUSE surveys (e.g. 



Rarhfoi'.l <>t al. 200*3) . In several low column 



density lines of sight ISavage et al.l (119771 ) did not quote an uncertainty for the H 2 columns. 
In these cases we use the uncertainty a{\ogN^ 2 ) = 0.20 as suggested by D. WeltjQ and i 
comparable to the maximum quoted value for all observations. 



is 



Many of the Nca co lumn s are from the Coper nicus ultraviolet absorption observations of 
Jenkins fc Shayal (119791 ) and IJenkins et al.l (Il983l). Additional C I columns are provided for 



example, by GHRS ob servations and FUSE (IZsargo et al.lll997l; IZsargo k, Federmanl 12003 
Jenkins fc Tripp! 120011 ). The method used by IJenkins &: Tripp! (120011 ) to determine C I 



columns does not yield a direct measure of the random errors in the observed velocity 
integrated column density. We adopt an unce rtainty comparable to simil ar STIS and FUSE 



observations at similar column densities (e.g.. ISonnentrucker et al.l 120021 ) . 



lit t p : / / astro . uchic ago . edu/ home/ web / welty / coldens . ht ml 
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MODELS 



3.1. Reaction Rates with PAHs and Carbon Chemistry 



The PAH rates i n IWolfire et al.l (120031 ) were calculated using the equations given by 
Draine fc Sutinl (119871 ) for a disk shaped PAH with the num ber of Carbon atoms Nq = 35 
and a PAH abundance of 6 x 10~ 7 per hydrogen nucleus. IWolfire et al.l (120031 ) found it 
necessary to modify the rates by a factor $pah to match a preliminary assessment of the 
observed N C i/N C u ratio in the diffuse ISM. The factor $pah is to include a wide range of 
unknowns in th e rates including the PAH size, geometry, and abundance. Although the full 
Draine fc Sutinl (119871 ) formalism is used in our code, over typical temperatures and densities 
found in the diffuse ISM the rates can be closely approximated by simple formulae. For C + 
recombination on PAHs we have 



C+ + PAH" -> C + PAH , jfei = 2.4 x IO^fah^ -0 ' 5 cm 3 s" 1 
where T 2 = T/(100 K), and for charge exchange 

C + + PAH -> C + PAH+, k 2 = 8.8 x 10" 9 $ PA h cm 3 s" 1 . 



(2) 



Several rates were given in IWolfire et al.l (120031 ) Appendix C2, mainly related to the charge 
balance for hydrogen and the electron abundance. We provide a summary of our rate set in 
Appendix [A] 

In addition to reactions with PAHs, C + can also recombine with electrons in the gas 
phase via radiative or dielectronic recombination. It has recently come to our attention that 
dielectronic recombination of C + can be significant in diffuse gas (G. Ferl a nd, p rivate com- 
munication). The total (dielectronic plus radiative) rates by lAltun et al.l (120041 ) and those 



posted on-line for lower temperatures by Badnell ( htt p://amdpp.phys.st r ath.ac .uk/tamoc/DATA/ ) 



are higher by a factor of ~ 2 at 100 K than rates by lNahar fc Pradhanl (Il997l ). For the ma- 
jority of this paper we use the Badnell rates, but discuss the implications if they are lower 
as suggested by the previous data. A fit to the total rate is given by 



C + + e -> C° k 3 = 1.8 x 10" n T 2 -°- 83 cm 3 s" 



(3) 



where the fit is good to within 8% for 20 K < T < 300 K. 

The C + is formed mainly by the photoionization of C by the interstellar radiation field 



C + hv -> C + + e, h = 2.1 x 10~ 10 G exp(-2.6A 



(4) 



where Go is th e FUV (6 eV < hv < 13.6 eV) interstellar radiation field mea sured in units 
of the iHabingl (119681 ) field (G = 1 is 1.3 x 10~ 4 erg cm -2 s -1 sr" 1 ). In IWolfire et al. 



- 6- 



( 120031 ). we adopted Go = 1.7, a field strength comparable to the iDraind (119781 ) interstellar 
field. Balancing formation and destruction processes, in clouds optically thin to the FUV 
radiation, the N C i/N C u ratio is mainly a function of n^p^/Go. 



3.2. H 2 Formation 



Molecular hydrogen formation in diffuse gas proceeds via reactions with atomic hydro- 
gen on grain surfaces. We adopt a formation rate per unit volume which goes as nn^iR where 
n is the hydrogen nucleus density, nni is the density of atomic hydrogen, and R is a constant 
of order 3 x 10~ 17 cm 3 s _1 . The constant includes the surface area of grains per hydrogen 
nucleus. (If the metallicity were to vary as Z, and the grain surface area per H proportional 
to Z ', then the rate R would go as Z.) Previous investigators have included both a gas and 
grain temperature dependence in this rate equation intended to include the collision rate of 
atomic hydrogen with grains, the sticking coefficient of hydrogen on grains (which is a func- 
tion of both the hydrogen temperature and grain temperature), and the thermal diffusion of 
atoms across grain surfaces and thermal evaporation from these surfaces (which is a function 
of the grain temperature). The collision rate increases with tempe rature while the stick- 
ing coefficient and formation efficiency decrease with temperature ( Hollen bach fc Salpeter 
197ll : iBurke fc Hollenbachl fl983l ; IPirronello eiTal~1ll999l : ICazaux fc Tielensl 120041) . In light of 
the uncertainties of the details of the various processes, iKaufman et al. fjl999h adopted an 
average rate which was independent of temperature. 



The dissociation rate of H 2 per unit volume goes as 

G In n2 p(N R2 ) exp(-2.5A v ) 



(5) 



where I is the unshielded photodissociation rate in the local interstell ar field (G n = 1), an d 
f3(N R2 ) is the H 2 self-shielding factor. We take I = 4.7 x 10~ n s" 1 from kbgrall et all Jl992h . 
The unshielded rate could change depending on the detailed population of the molecular 
hydrogen rotational and vibrational levels resulting from different incident FUV fields, gas 
densities, and temperatures. We have compared our constant unshielded dissociation rate 
I, with res ults produced with t he code available on-line from the Meudon GroupB As 
discussed in lLe Petit et al.l ( 120061 ) they include a more detailed treatment of H 2 dissociation 
and shielding. We ran their code at n = 30 and 100 cm -3 , Go = 1.7, and Ay = 0.1, 0.3, 
and 1, and found good agre ement with less than a 10 % difference in values of I. We use 
the self-shielding factor from iDraine fc Bertoldil (119961 ) (their equation 37) which is a fit to 



http://aristote.obspm.fr/MISI 
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their detailed H2 destruction calculation. The equilibrium molecular fraction 2nH 2 /(^Hi + 
2nH 2 ), for unshielded H 2 (/? = 1) equals 2nR/G I. For constant /, the molecular fraction 
is proportional to 2hR/Gq. Thus, to first order and ignoring shielding by dust or H 2 , the 
Nqi/Nqu ratio depends on ^p\nn/G while the N^/Nm ratio depends on Rn/G . 



3.3. 2-Sided Models 



We see from sections 13.11 and 13.21 that the Nqi/Nqh ratio and the H2 column A^h 2 
depends on $pah, R, and the ratio of Go/n. The total column density iV = Nm + 2Nn 2 
also enters due to the extinction of the incident FUV field in equations (TjJ and fl5J) and also 
in the self-shielding factor f3 (which depends on the H 2 column density). The total column 
density is known for each source. Thus we can construct models for each source by holding 
the total column density fixed, and varying Go/n, $pah, and R, until good fits are obtained 
to the observed columns densities, Nqi/Nqh and 2N^ 2 /N 

Since many of the observations have cloud columns Ay ^ 1 the effects of radiation in- 
cident on the near and far side of the cloud need to be included. This radiation affects both 
the photoionization of C and the dis sociation of H 2 . We have modified the 1-sided Pho- 
todissociation Region (PDR) code of Kaufman et al.l (119991 ) to include the 2-sided incident 
radiation. The code calculates the equilibrium chemical abundances and gas temperature in 
a gas layer exposed to X-rays and FUV radiation. In the one-sided case, calculation of the 
chemistry, cooling, and line transfer proceeds from the surface to the cloud center in a single 
pass since these parameters depend only on the cloud properties closer to the surface. In the 
two-sided case, at a given point in the cloud, the optical depth in the cooling lines and the 
H 2 column is required towards both surfaces, thus an iterative procedure is required. We 
first calculate the structure of a cloud illuminated from one side to a depth of one-half the 
visual extinction of the desired two-sided model. This provides us with initial estimates of 
line optical depths as well as shielding columns of H 2 and CO. We then iterate the chemical 
abundance profiles, optical depths, and shielding factors. 



We have further modifie d the physics and ch emistry of the iKaufman et al.l (119991 ) code 
according to the discussion in lWolfire et al.l (120031 ) who assumed a higher abundance of PAHs 
and a slower rate of interaction of ions and PAHs. These changes largely offset each other, 
with a resulting minor effe ct on the g r ain h eating rates and ion chemistry. In addition, we 
have adopted the results of Peauignotl (119901) fo r the rate coefficient for collisional excitation 
of O I by H, and those of iMcCall et al.l (120031 ) for the H^" dissociative recombination rate 
coefficient. We have also incl uded the illum i natio n of the cloud by the interstellar soft X-ray 
radiation field as discussed in lWolfire et al.l (120031 ). Results from this 2-sided code have been 
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previously presented in lNeufeld et al.l (120051 ). and lSnow fc McCalll (120061 ). 



We have compared our calculated H 2 column densities with those produced by the 
Meudon code. For the same model parameters in section 13.21 and running the code in 2- 
sided mode, we find that our H2 columns agree to within a factor of 2 (with our columns 
lower) and with the largest difference at low column densities (Ay ~ 0.1). The agreement is 
within the typical observational error. 



3.4. Fitting $pah and R 

We determine the best fit $pah and R by computing a grid of models consisting of 12 
values of $pah (0.01, 0.1, 0.2 ,0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1) and 11 values of R (0.5, 
1, 2, 3, 4, 5, 6, 7, 8, 9, 10) xl0~ 17 cm 3 s _1 . At each grid point, (i.e., for fixed $ PA h and R), 
we find the best fit Go/n values for each source by minimizing the xl 

log /£ odcl - log / H2 1 2 , / log f^ dcl - log fc I 



I ^[log/Ha] J I O-Pog/ci] J 

where xl is the x 2 value for source s, /h 2 and fci are the observed values for the column 
density ratios fu 2 = 2N R J[N U1 + 2Nn 2 ], fci = N C i/N C i h /f 2 odel , and f$? del are the calcu- 
lated values of /h 2 and fci (dependent upon N - which is measured, and on the unknowns 
$pah, R, and G /n). 

We will calculate the x 2 value over a subset of the sources within a range of Ay val- 
ues. Thus the reduced x 2 value at grid point ($pah, R) is the sum over sources within the 
restricted Ay bin: 

X 2 ($PAH^)=$> 2 /(iV s -2). (7) 

s 

where N s is the total number of sources in the bin. 

We find the uncertainty in the measurements cr(log/H 2 ) and cr(log/ci) using the quoted 
uncertainties <r(log Nm), cr(logiVH 2 ), cr(logA^ci), for the H I, Ho , and C I columns, along 
with the standard propagation of error formula (e.g.. lTaylorlll997l ). Thus 

a(log/ H2 ) = K 1 °g Ar H 2 ) 2 + cr(logiV) 2 ] 1/2 , (8) 
where cx(logiV) is the uncertainty in the total column density, 

a(logiV) = [a(logiVH 2 ) 2 ^ 2 +cr(logiVHi) 2 ^i] 1/2 , (9) 
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and /hi is the fraction of atomic hydrogen /hi = Nm/[Nm + 2Nu 2 ]. For observed C II 
column densities, the uncertainty in log/ci is given by 

t7(log/ci) = [a(logiVci) 2 + a(logiV CII ) 2 ] 1/2 , (10) 

while for estimated C II column densities the uncertainty in log/ci is given by 

a(log/ CI ) = [<7(logiV CI ) 2 + <7(logiV) 2 ] 1/2 . (11) 

In practice, for the x 2 test, we set Go = 1.7 and vary n. We estimate the initial range 
of n values from the analytic expressions in Appendix [B] We run models stepping in values 
of n until we cross the minimum in \ 2 - Then we apply a parabolic interpolation to find the 
minimum \ 2 and associated values of density and column densities. 



4. RESULTS and DISCUSSION 

4.1. Fixed Column Density Models 

We first carry out model runs for four fixed cloud column densities Ay = 0.1, 0.5, 1.0, 
and 1.5 with variable $pah, R, and Gq/u. The data and model results are displayed in plots 
of Nqi/Ncii versus /h 2 with source column densities grouped in four Ay bins. These plots 
are intended to show the range of observed values, the effects of varying $pah and R, and 
an "eye ball" assessment of reasonable values. The Ay bins span the range of observed Ay 
values with the smallest being Ay = 0.03 and the largest Ay = 2.13. The Ay bins have been 
chosen to display natural groupings of the data which we further test using the \ 2 results. 

We present in Figure [1] the observed and calculated ratio of Nci/Nqu versus /h 2 for 
observations in the range 0.03 < Ay < 0.25 and N U2 < 10 17 cm 2 . Figure [1] focuses on 
the low /h 2 ( ~ 10~ 4 ) and low Ay directions. The curves are for a model cloud column 
density of Ay = 0.1, with Gq = 1.7 and 5.1, and density points at n — 10, 20, 30, and 
40 cm -3 (with the exception of two curves with three density points at n = 5, 10, and 20 
cm -3 ). Diamonds (O) indicate models with Go = 1-7 and triangles (A) indicate models with 
G = 5.1. We have labeled the model curves with their corresponding values of $pah/-R-i7 
where i?_i 7 = .R/lxlO -17 . The curves with $ PAH /#-i7 < 1/3 use R = 3x 10 -17 cm 3 s -1 with 
values of $ PA h = 1 ($pah/-R-i7 = 1/3), $pah = 0.5 ($ PA h/-R-i7 = 1/6), and $ PA h = 0.01 
($pah/-R-i7 = 1/300). We have added three additional curves with lower values of R namely 
R = 2x 10- 17 cm 3 s" 1 , $p AH = 1 ($pah/#-it = 1/2); R=lx 10~ 17 cm 3 s" 1 , $ PAH = 0.5 
($pahAR-i7 = 1/2); and R = 1 x 10~ 17 cm 3 s _1 , $pah = 1.0 ($pah/^-i7 = 1) and one 
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additional curve with higher values of $pah namely R = 2 x 10 17 cm 3 s l , $pah = 2.0 
($pahAR-i7 = !)• 

The curves with the same value of $pah/-R-i7 ; although not exactly the same, differ by 
less than a factor of 1.5. We see that values of the ratio $pah/-R-i7 ~ 1/2 — 1/3 are required 
to match the bulk of the observations. For $pah < 1 this requires R < 2 x 10 -17 cm 3 s _1 
and an H2 formation rate smaller than the standard 3 x 10~ 17 cm 3 s™ 1 is necessary. 

The subsequent Figures (TSJ El and 0J focus on the high molecular abundance lines 
of sight (/h 2 > 10~ 3 and Nn 2 > 10 18 cm -2 ). The observations are limited to the range 
0.25 ^ Ay < 0.75 (Fig. 0), 0.75 ^ Ay < 1.25 (Fig. [3]), and 1.25 ^ Ay < 2.13 (Fig. @]) 
and model curves are for Ay = 0.5 (Fig. [2j), Ay = 1.0 (Fig. [3]), and Ay = 1.5 (Fig. [4]). 
The models all use R = 3 x 10~ 17 cm 3 s _1 with values of $pah = 0.01, 0.5, and 1.0, and 
density, n = 10, 30, 100, and 300 cm -3 . Comparing Figures [T] and [2] we see a large jump in 
the molecular fraction with no observations at intermediate values. This jump is probably 
the result of the abrupt turn on of molecular hydrogen due to the eff ects of self-shielding . 



Self-shielding of H 2 at iV > 5 x 10 20 cm" 2 (Ay > 0.25) was seen in thelSavage et all (119771 ) 



data and also in the more recent FUSE surveys (jGillmon fc Shuilll2006l ). In Figure [H two 
obvious outliers are 23 Ori (Ay = 0.28) at high Nqi/Ncii ratio (= 1.1 x 10~ 2 ) and rr Sco 
(Ay = 0.28) at low iV C i/iVcii ratio (= 1.2 x 10~ 3 ). These are excluded from the \ 2 tests 
and discussed further in subsection 14.51 

The variation of the curves shown in Figures [T] through H] as functions of Go/n, 0pah 5 
and R is discussed in Appendix [B] Points with the same Go/(nR) and fixed column, or 
Ay, have the same molecular fraction /h 2 independent of $pah- At the highest values of 
Go/n the molecular fraction is lowest due to rapid photodissociation of molecular hydrogen. 
In the limit of no H 2 self-shielding and no dust extinction, the limiting molecular fraction 
is given by fu 2 = 2nR/(G I). As G /n decreases, the molecular fraction increases with a 
limiting value of /h 2 — 1 f° r fully molecular gas. As the cloud column density (and thus 
Ay) increases the local FUV field drops due to dust extinction and the molecular fraction 
rises due to a lower H 2 photodissociation rate. In addition, the self-shielding by H 2 raises 
the molecular fraction as the column density increases. 

The variation in the Nqi/Nqii ratio can be understood in four limiting regimes: high/low 
$pah, and high/low Go/n. First consider the high $pah (0.5 and 1.0) curves. For low 
Go/n the destruction of C + is dominated by recombination on PAH~ and Nqi/Nqu oc 
$p AH (n/G ) 2 (See Appendix [C] eq. |Clj ). For high G /n, the destruction of C + is dominated 
by charge exchange with PAH and electron recombination and thus Nqi/Nqu oc ^pah^/Go 
(eq. jC2j ). In the case of low $pah = 0.01 the C + destruction is dominated by recombination 
with free electrons n e . For low Go/n the electrons are provided by photoionization of C, and 
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Nci/Ncn oc n/Go (eq. jC3j ). For high Gq/u, i.e., low values of n, the electrons are provided 
by the soft-X ray ionization of H, and Nqi/Nqu ^ /Go (eq. |C4j ). In the molecular 
hydrogen case, the molecular fraction increases with column density due to the drop in FUV 
field and larger H 2 self-shielding. However, in the case of the Nqi/Nqu ratio, it is just the 
dust extinction that leads to the increase in the ratio with column. 



4.2. x Results 

We next present the \ 2 results obtained by binning the sources in four Ay bins. But 
note that for each source model, we fix Ay to the observed column. Figure shows the 
results for sources in the range 0.03 £ A Y < 0.25 and Nu 2 < 10 17 cm . There appears to 
be two local minima: one at low R values with R = 1 x 10~ 17 cm 3 s _1 and 0.2 ^ $pah ~ 0.4 
and one at R ~ 2 x 10~ 17 cm 3 s^ 1 and <3?pah ~ 1- The Nqi/Nqu ratios in this Ay range 
are ^ 10~ 3 and both equations (1B15P and (IB 17ft apply. At low $pah the C + destruction 
by electron recombinations dominate (right hand side of eq. IB17P and the dependence on 
$pah drops out. Moving to the upper right in the \ 2 plot corresponds to lower densities and 
higher temperatures. There, electron recombination is less effective in producing C I, the C I 
abundance drops, and the fits become poorer. At sufficiently high $pah 5 the PAH reactions 
dominate, the C I abundance rises, and the fits improve. We shall demonstrate that the low 
R and $pah minimum corresponds to gas densities that are closer to those expected for the 
diffuse ISM and are therefore favored. The minimum Xmm = 1-6 and lies within a closed set 
of contours. 

We next show the results for sources in the remaining Ay ranges 0.25 ^ Ay ^ 0.75 
(Fig. 10 sources), 0.75 < Ay < 1.25 (Fig. 0: 12 sources), and 0.25 < Ay < 2.13 (Fig. 
[8]: 28 sources). In these \ 2 plots we do not consider the exact position of the minimum to 
be as significant as in Figure [5] since it generally lies in a trough of low values within open 
contours. An approximate fit to the minimum trough is given by 

$ PAH = 0.22 ( — — — ) - 0.40 R > 3 x 10~ 17 cm 3 s" 1 . (12) 



10 17 cm 3 s 



and shown in Figure [H 



The x 2 contours shows that R and $pah are correlated. The minimum lies at roughly 
a 45 degree angle with $pah proportional to R. For a given observation, that is a fixed /h 2 , 
and Ay, from equation (1B9I) . the ratio Ru/Gq must be constant along the line of minimum 
X 2 . In addition, for a fixed Nqi/Nqu, from equations (IB15I) and flB17j) the ratio Go/ (ti^pak) 



must be constant along the line of minimum % 2 . Together these imply that i?/$pAH must 
be constant or R oc ^pah- 
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There is a suggestion from our figures that the best fit $pah and R increases with column 
density. For 0.25 ^ Ay ^5 0.75 (Fig. [6]) the minimum x 2 hes at R = 3 x 10 -17 cm 3 s _1 , and 
^pah = 0.2, while for 0.75 ^ Ay ^ 1.25 (Fig. [7j), there is a minimum trough that extends 
from R = Ax 10~ 17 cm 3 s _1 , and $ PAH = 0.4 to R = 6 x 10~ 17 cm 3 s _1 , and $ PA h = 1-0. If 
this effect is real then it is not entirely clear what is the physical reason behind it. Perhaps 
shock processing along low Ay lines of sight tend to reduce grain surface area and produce 
smaller R (and therefore smaller <3>pah since $pah oc R). The statistics however, are rather 
poor for separate Ay bins and additional observations are required to confirm the trend. 

Note that as the H2 formation rate R increases the density n must decrease (for fixed 
Go) in order to maintain the same molecular fraction. At some large value of R the densities 
become much lower than typically observed in the diffuse ISM. In Figures [5] and [H] each R 
and $pah grid point has several sources contributing to the \ 2 value and there is a best fit 
density for each source. We show in Figures [9] and [10] the average density (n), obtained by 
averaging the best fit density over the sources contributing to each R and $pah grid point. 
Recall that in computing the x 2 grid we are finding the best fit ratio Go/n and have fixed 
Go at 1.7. Thus the density could be higher for higher Go- Also shown in Figure [10] is 
the line along the x 2 minimum trough (eq. [12]). The density plot (Fig. [TUT) combined with 
the x 2 pl°t (Fig- E]) shows that the typical diffuse cloud (n ~ 30 cm -3 , G ~ 1.7) lies at 
R ps 3.5 x 10~ 17 cm 3 s _1 and $pah = 0.4. For the low column densities, the contours centered 
on R ps 1 x 10~ 17 cm 3 s -1 , and $pah = 0.2 — 0.4 correspond to gas densities closer to that 
expected in diffuse gas and we therefore favor these over the larger R and $pah values. 



4.3. Effects of C + Recombination Rate 



The recombination rate of C + with electrons is most important at low 0pah or high 
values of Gp/n (se e Appe ndix IB]) . The total rates (including dielectronic reco mbination) of 
Nahar fc Pradhan] (119971 ) are higher by a factor of ~ 2 at 100 K compared to lAltun et al. 
( 120041 ) and further diverge at lower temperatures. The Nqi/Nqh ratio can thus vary by 
a factor of ~ 2 when C + recombination with electrons dominates the production of C I. 
In previous sections of this paper we have used the higher rates, but since the rates are 
uncertain, we have also investigated the effects of the lower rates. Figure [11] shows the 
resulting x 2 plot for the range of column d ensities 0.25 < Ay < 2.13 and A^ H2 > 10 18 



cm 



The effects of the iNahar fc Pradhan! (119971 ) rates more narrowly constrain the best values of 
0pah and R by excluding the parameter space of low 0pah and large R. This is because of the 
decreased importance of the electron recombination rates and thus the observed Nqi/Nqii 
ratio can only be achieved through higher PAH reaction rates. The best fit minimum trough, 
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however, is not changed significantly compared to Figure with only a slight extension to 
lower R ~ 2.0 x 1(T 17 cm 3 s" 1 and PAH ~ 0.2. 



4.4. Ay < 0.25 Lines of Sight 

For our low column density lines of sight, the combination of a low molecular hydrogen 
fraction, / Ha and a relatively high N C i/N C n ratio indicate that R must be low or at least 
a low ratio of -R/$pah (ref. Figs. pQ and [5]). Are there unusual physical conditions, or 
additional destruction processes that are acting in these regions which could account for 
the low H 2 columns? We first consider the effects of warm neutral medium (WNM) along 
the line of sight. The WNM consists of low density (n ~ 0.3 cm~ 3 ) warm (T ~ 8000 K) 
gas. Such gas will likely contain C II and H I, but very little C I and H2. This is because 
the formation rates of C I and H2 are proportional to n 2 , while the destruction rates are 
proportional to n and thus the abundance drops in lower density gas. In addition, at high 
temperatures the C + recombination rates fall. The H 2 formation rate may also be lower in 
the WNM because at high temperatures the H I atoms colliding with grains will more likely 
bounce than stick ( Burke fc Hollenbachlll983l ). We tested results for an Ay = 0.1 cloud with 



Go = 5.1 and n = 30 cm -3 in which half the column is WNM gas. Here we have included 
the density and temperature effects in the rates but not reduced R due to a smaller sticking 
coefficient (see Fig. [121 arrow W'). The model point which includes WNM moves towards 
lower Nci/Nqii ratio (by a factor of ~ 2) and lower /h 2 (by a factor of ~ 10). Part of the 
drop in /h 2 is due to lower H 2 self- shielding in the n = 30 cm -3 cloud. The resulting vector 
moves away from the low Ay observations, and along a path which is nearly parallel to the 
standard models. Because the addition of WNM along the line-of-sight only moves points 
parallel to the standard models, it is unlikely that the presence of a WNM component can 
explain the /h 2 an d Nqi/Nqh results at low column. 

We next consider the effects of a higher cosmic ray ionization rate along low column 
density lines of sight. The effects of increasing the primary ionization rate by a factor of 10 
has little affect on the H2 columns since photodissociation dominates the destruction of H2. 
The iVci/iVcii ratio increases slightly (by ~ 30%) due to the greater electron abundance from 
cosmic ray ionizations (Fig. [12] arrow 'C'). We conclude that neither WNM nor enhanced 
cosmic ray rates can explain the low Ay observations. 

Another possibility is that the low Ay lines of sight are strongly affected by grain 



processing in interstellar shocks. iJones et al.l (119961 ) demonstrated that sputtering and grain 



grain collisions can vaporize grains, and reduce the grain surface area in sufficiently fast 
v s 150 km s -1 shocks. In the interstellar medium, shock speeds tend to scale with n -1 / 2 , 
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so that lower density gas will experience higher velocity shocks. Low Ay gas may have higher 
fractions of low density WNM gas. We note however, that shocks of speeds v s ~ 100 km 
s -1 tend to shatter the larger grains and produce a large population of small grains which 
increases the surface area. We speculate that another type of grain processing occurring in 
shocks is to sputter the surfaces "clean" and thereby modify the characteristics of the atomic 
hydrogen adsorption sites (e.g., binding energy). 

A final possibility is that the low Ay columns are lines of sight which graze the edges of 
larger clouds. In this case, the typical shielding column density (which is weighted towards 
the smallest column direction, i.e., perpendicular to the line-of-sight) is much less than the 
line-of-sight column that we measure. We examine this possibility with a simple test. We 
take a cloud with a total Ay = 1 through the center and consider a ray which passes through 
the cloud with a column of Ay = 0.1. At the midpoint, such a ray passes through a minimum 
depth of Ay = 2.5 x 10~ 3 measured from the cloud surface. We take the nn 2 , nni, nci, and 
no ii abundances as a function of cloud radius from model runs for Ay = 1, Go = 5.1, and 
n = 30 cm -3 . We then integrate the Nm, Nu 2 , Nci, and Nqii column densities along the 
grazing line of sight. Results are shown in Figure [12] as an arrow labeled 'G'. We see that 
compared to an Ay = 0.1 cloud, the /h 2 decreases due to the diminished H 2 self-shielding at 
the cloud edge. Only a small molecular column exists between the FUV illuminated surface 
and points along the ray and thus the self-shielding is ineffective and the H2 is dissociated. 
We also see that the Nci/Ncn ratio rises. This is because the FUV field is lower by a 
factor of two compared to the Ay = 0.1 cloud since the radiation can not penetrate from 
the backside of the larger cloud. The lower field reduces the photoionization rate of C to C + 
and increases the abundance of PAH - which converts C + to C. (The lower field also reduces 
the H 2 dissociation rate but this reduction is offset by the lack of H 2 self-shielding). The 
results significantly move the model points towards the observations. We consider grazing 
incident lines of sight a potential explanation for the low column density observations without 
resorting to unusual grain properties. 

For Ay ^5 0.25 we find two lines of sight to be anomalous. These are e Per (Ay = 0.16, 
/h 2 = 0.20) and 59 Cyg (Ay = 0.11, /h 2 = 0.19). These lines of sight have molecular 
fractions which are ^ 10 3 times greater than others with similar (low) column densities and 
require different rates or environments than the low /h 2 cases. Including e Per and 59 Cyg 
in the x 2 tests for the 0.03 ~ Ay ^ 0.25, Nu 2 < 10 17 sources raises the minimum \ 2 from 
1.6 to 9.2 and thus these sources can well be considered to be part of a separate population 
than the rest of the low Ay and low /h 2 lines of sight. The k Ori line of sight at slightly 
higher column density (Ay = 0.17) is part of the low f H2 group as are e Ori (Ay = 0.14) and 
15 Mon (Ay = 0.13). One explanation for these peculiar sight lines could be that shocks 
of speeds v s ~ 100 km s -1 have significantly increased the population of small grains. The 
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resulting increase in grain surface area will increase the FUV extinction (and decrease the H2 
photodissociation rate) and also increase the H 2 formation rate. We find that with the FUV 
extinction and H 2 formation rate a factor of three higher than normal we can reproduce the 
Nci/Nqu ratio and molecular fraction for the e Per line of sight. Unfortunately, there are 
no low resolution IUE data for e Per to examine the FUV rise in the extinction curve. We 
also note that 59 Cyg is a variable Be star and column densities for this line of sight are 
highly uncertain (D. Massa, private communication). 



4.5. A v > 0.25 Lines of Sight 



At columns densities Ay > 0.25 we found two points to be anomalous: 23 Ori and it Sco. 
Weingartner fc Draind (120011 ) modeled the line of sight towards 23 Ori and found that they 
had some difficulty matching the high Nci column even when recombination on grains was 
included. They proposed several additional processes which might be working to enhance 
the C abundance including the dissociative recombination of CH + . The CH + abundance can 
be enhanc ed compared to normal equilibrium chemistry due to the effects of non-therma l 
chemistry (IFlower fc Pineau des Foretslll998l ; IJoulain et al.lll998l ; IZsargo fc Federmanll2003l ). 
If turbulent velocities are sufficiently large (v ^ 3 — 4 km s _1 ) and turbulent dissipation can 
drive reactions with temperature barriers, then a large column of CH + can be produced 
djoulain et aDll998h . 



Another process which might account for the anomalous columns is the time dependence 
of the chemistry. The line of sight towards 23 Ori shows a large Nqi/Nqh ratio at small 
molecular fractions fn 2 compared to the equilibrium abundances. Consider a parcel of gas 
which has been recently shocked (i.e., H 2 dissociated and C I ionized). The molecular 
hydrogen is produced at a rate of ~ 3 x 10~ 17 n 2 s _1 . The time to double the current fractional 



abundance nu 2 /n ~5x 10~ 3 is about t ~ 5 x 10~ 3 /(3 x 10~ 17 n) s or t ~ 5.3 x 10 4 yr for 
n = 100 cm -3 . On the other hand, the time to double the current ratio Nqi/Nqh ~ 0.01 is 
about t ~ 0.01/($pah2.4 x 10~ 7 n P AH-) s or t ~ 62 yr for n = 100 cm -3 , where we have used 
C + recombination with PAH - , $pah = 0.5, and npAH- from equation (1B14|) . The time to 
reach the equilibrium Nqi/Nqh ratio is much shorter than the time to reach the equilibrium 
/h 2 fraction. Thus, for a recent (t ^5 5 x 10 4 yr) passage of a shock, the Nqi/Nqii ratio has 
had time to reach equilibrium while the H 2 fraction is dissociated compared to equilibrium. 

For 7r Sco consider a parcel of molecular gas which has been recently illuminated by 
a source of FUV radiation. The time to dissociate the current /h 2 ~ 0.01 at cloud center 
is t ~ (GqIP[Nb. 2 ] exp[— 2.5Av])~ 1 or t ~ 1.1 x 10 6 yr, where we have used G = 5.1, 
N R2 = 1 x 10 19 cm" 2 , (3[l x 10 19 cm- 2 ] = 1.7 x 10~ 4 (Eq.[B3]), and A v = 0.13. The time to 
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ionize the current N C i/N C n ~ 10~ 4 ratio is about t ~ (2.1 x 10~ 10 G exp[— 2.6v4 v ]) _1 s or 
t ~ 41 yr. Thus for a molecular clump illuminated within the past ~ 10 6 yr, the Nqi/Nch 
ratio is in equilibrium while the H 2 fraction is much higher than its equilibrium value. 

One complicating factor which we have not yet considered is the effects of multiple 
clump components along the line of sight. In general, taking the observed column density 
and splitting it up into multiple clouds exposes more surfaces to the interstellar radiation 
field. In the case of optically thick clumps, this will decrease the N C i/N C u ratio and decrease 
the molecular fraction / Ha . For the case of optically thin clumps, the Nqi/Nqh ratio will 
not change but /h 2 will drop due to reduced self-shielding. However, the effects of multiple 
clumps is important only for components of approximately equal column density. If one 
component dominates then the remaining smaller components do not contribute much to the 
total column. Since in general, the UV absorption spectroscopy observations are dominated 
by a single component, our neglect of multiple clouds should not significantly alter our 
results. 



5. Summary 

We have carried out an analysis of observations of the C I, C II, H I, and H 2 column 
densities in the diffuse ISM towards 42 lines of sight spanning a range of column densities 

0. 03 < Ay < 2.13. We have fitted these columns using a 2-sided PDR model which simul- 
taneously solves for the thermal and chemical balance throughout the layer. The relative 
columns of H I and H 2 are dependent on the rate coefficient R for H 2 formation on grain 
surfaces, on the ratio Go/n, and on the total column N (or Ay) of the illuminated cloud. 
The relative columns of C I and C II depend on the electron and PAH abundances, the rate 
coefficients for neutralization of C + with PAHs (or $pah), the total column in the illumi- 
nated cloud, and Go/n. Observations set the total columns of each line of sight plus the 
fractional abundances of H, H 2 , C, and C + . Our PDR models self consistently determine 
the electron, PAH + , PAH, and PAH - abundances for each model. We have performed a \ 2 
test using R, $pah, and Go/n as our free parameters to find the best fits to the fractional 
abundances along each line of sight. 

The results are shown in Figures 151 through [T0| where we show the \ 2 contours in the R, 
$pah plane, given the best fitting G /n for each individual line of sight. Note that Figures 
151 through [TUl are the x 2 f° r these values of R, $pah for all 42 sources. Our main conclusions 
are as follows: 

1. For low column (Ay < 0.25) lines of sight, and for the low molecular columns and fractions 
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that accompany them (Nu 2 < 10 17 cm 2 , /# 2 



< 10" 4 ), we find R ~ 1 x 10~ 17 



cm 3 s 



PAH 



0.2 - 0.4. 



2. At higher column densities (0.25 < Ay < 2.13), we find R 
$pah ~ 0.4. 



3.5 x 10 



-17 S 

cm s 



and 



and 



3. Our x 2 fits show a correlated range of $pah and R values that provide as good fits as our 
"best fits" quoted above. For A v > 0.25, we find good fits for ® PAH = 0.22(i?/10" 17 cm 3 s" 1 )- 
0.40 as long as R > 3 x 10~ 17 cm 3 s _1 . Very high values of R (or therefore $pah), however, 
can only be fit with very low average densities in the clouds along the 42 sight lines. Figure 
[TUI shows that for densities (equivalently thermal pressures) thought to be typical of local 
diffuse interstellar clouds, a value of R ~ 3.5 x 10~ 17 cm 3 s _1 (and therefore $pah ~ 0.4) 
is preferred. This value is consistent with previo us determin ations of R made by comparing 
observations of H 2 with model calculations (e.g. 



Jural 11975( 1 . 



4. The low values of R for low column density sight lines may indicate shock processing 
that decreases grain surface area or modifies grain surfaces along these sightlines, thereby 
reducing the H 2 formation rates, or a line of sight which grazes a larger cloud. 

5. Our resu l ts for $pah indicate that the crude estimates of these ratios adopted by 



Wolfire et all (120031 ) were roughly correct. They confirm the importance of PAHs in de- 



termining the ionization level of atomic ions (in this case C vs C + ) in interstellar clouds. 

M.G.W. was supported in part by a NASA Long Term Space Astrophysics (LTSA) grant 
NNG05GD64G. We thank an anonymous referee for a careful reading of our paper and for 
comments that improved the presentation. 



A. Chemical Rates 



The chemical rates with PAHs were presented in IWolfire et all (120031 ) Appendix C2 
for reactions involving the ionization balance of hydrogen. Here we reproduce the rates 
along with the Carbo n rates used in this paper. Our starting point are the rates from 
Draine &: Sutinl (119871 ) with N c = 35 carbon atoms, and disk PAHs with a radius a = 



(N c / 1-222) a5 . We use a PAH abundance of n PAH = 6 x 10~ 7 n. IWolfire et all ((20031) intro- 



duced the factor $pah to account for a wide range of unknowns in the PAH distribution and 
rate coefficients. Simple fits to the rates appropriate to diffuse ISM conditions are presented 
in Table [1 
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B. Estimating Radiation Field Go and Density n 



In this appendix we derive simple expressions to estimate the FUV radiation field Go 
incident on the cloud and the cloud density n from the observed column densities of iV(HI), 
-Wh 2 > N C i, and N C n. In the limit of no dust extinction, the equilibrium abundance of H 2 is 
given by 

G /% 2 /3[iVH 2 ]=i?nnHi , (Bl) 

where G is the incident FUV radiation field, I the unshielded photodissociation rate, P[N^ 2 ] 
is the H2 self-shielding factor, and R is the formation rate on grains. For H2 columns ~ 10 14 
cm -3 , /3[iVH 2 ] = 1 and the molecular fraction is given by 



H 2 



2nR 



(B2) 



To estimate Go and n we use a simple expression for the H2 self-shielding formula from 
Draine fc Bertokhl Jl996h 

N 



3/4 



where N = 1 x 10 14 cm -2 is the H 2 column where line self-shielding starts to become 
significant. We can integrate equation fIBll) through the cloud 



G IN> 



3/4 



dN, 



H 2 



Wo + iV H J3/4 

Letting x = N + Nn 2 the integral becomes 



Rn / dN m 



G IN 



3/4 f"o+Nn 2 dx > 



,'3/4 



'JV x 

Carrying out the integral and solving for Gq/u we have 

G 0.25iV H i 

n 



RnNn] 



R 



[N + AT H2 ]i/4 - iV 1/4 } N 3 /4 1 



or 



with 



G v R 
n I 



0.25N H] 



[N + N U2 }^ - N 1/4 \ N 3 /4 



(B4) 



(B5) 



(B6) 

(B7) 
(B8) 
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For Nr 2 > N from equations (IB 71) . (1B8|) and lHollenbach fc Tielena (119991 ) the molecular 
fraction is given by 



Rn \ 4 / V x 



^ 2 UvJ UJ • (B9) 

For clouds where the optical depth of the dust becomes significant the local FUV field within 
the cloud is less than the incident field. We assume the flux on each side of a slab is 1/2Gq 
and the mean field, (Go) within the slab is approximately at a depth of 1/4 the total column 
density or at a depth of 1/4 Ay. Thus the mean field is given by 

(Go) = \g^-'^ + ±G e<-Wv) , (BIO) 

where a = 2.5 converts the FUV extinction at visual wavelengths to optical depth at the 
FUV dissociation energies of H2. Solving for Go we have 



2.5 

4 



1 + e ~ Av 



Setting G in equation flB7j) to the mean field (G ) we have 



. (B12) 

n 14- e ~~ Av 1 

The analytic solution is approximate since l/4Ay is an approximation to the depth 
where the local field equals the mean field. Comparing equation ( 1B12I) with our numerical 
results we find we can achieve better agreement with an additional factor of 1/2 on the right 
hand side. 

A second estimate of Gq/ti can be obtained from the Nqi/Nqh ratio. There are several 
destruction processes for C + to consider depending on the PAH~ abundance. If destruction 
occurs primarily through recombination on PAH" then the equilibrium abundance of C° is 
found by balancing recombination with photoionization of C 

G Icn C o = ^PAH^inc+npAH- > ( B13 ) 

where GqIc = Gq2.1 x 1CT 10 s _1 is the photoionization rate of C°, and $pah^i = $pah 2.4 x 
10~ 7 T 2 -0 ' 5 cm 3 s _1 is the C + recombination rate with PAH - . From I Wolfire et all ( 20031 ). the 
PAH - abundance in the regime ^pAH-/ n PAH° < 1 is given by 

rapAH- = 4.8 x 10 -5 n e n$pAHG " 1 cm" 3 , (B14) 
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with the electron density n e and hydrogen nucleus density n in units of cm 3 . Substituting 
this expression for tt-pah- into equation flB13j) . and solving for G /n we have 



^=4.1xlO- 3 $ PAH T : r a25 



n 



ci 



en 



-1/2 

(B15) 



where we have used an electron abundance of n e ~ 3 x 10 _4 n cm -3 (including metals, 
especially ionization of C, and ionization of H by soft X-rays), and substituted the observed 
ratio Nci/Ncn for n c o/n c +. 

If the abundance of PAH - is low enough, then the destruction of C + occurs through 
charge exchange with PAH or through electron recombination in the gas phase. The equi- 
librium abundance is then given by 

G Q Icn c o = nc+{k 2 npAno + k 3 n e ) (B16) 

with k 2 = 8.8 x 10 -9 $pah cm 3 s -1 , and k 3 = 1.8 x 10 -11 T 2 -a83 cm 3 s -1 . Using the PAH 
abundance of ~ 0.7 times the total PAH abundance (npAH — 0.7 x 6 x 10~ 7 n; Wolfire et al. 
2003), and an electron abundance of n e « 3 x 10~ 4 n cm -3 , we can solve for G /n 



— = (1.8 x 10- 5 $ PAH + 2.6 x 10- 5 T 2 -°- 83 ) 
n 



N, 



ci 



.iVcn. 

We find that the destruction of C + by PAH U dominates at a limit of approximately 

r > J-2 V PAH 



(B17) 



Go > T 2 1/2 ^ 
n ~ 1.2$p AH + 1.6T 2 _ 



0.83 



(B18) 



From equation ( 1B15I) . and setting T 2 = 1, and $pah = 0.5, this limit corresponds to an 
observed column density ratio Nqi/ Nc ii < 3 x 10~ 4 . Thus for N C i/N C n ~3x 10" 4 we use 
equation flB17j) while for N C i/N cu >> 3 x 10 -4 we use equation flB15j) . 



Similar to the Go/n estimate based on H 2 columns we also correct the incident FUV 
field for dust extinction within the cloud by setting the Go i n equations (IB15I) and (IB 171) 
to the mean FUV field (Go) an d solve for the extinction corrected Go/n. In solving for 
the minimum % 2 , we use the estimates for Go/n presented here along with a fixed value of 
Go = 1.7 to provide a first guess range of n. If the minimum in \ 2 is n °t crossed within the 
first range of n we take a second range given by n max x (n max /n min ) if higher densities are 
required or n min x (n max /n min ) _1 if lower densities are required (where n max and n min refer to 
the maximum and minimum densities in the range). The minimum in x 2 is generally found 
within the first or second range. We then do a quadratic fit of the \ 2 values as a function of 
density to find the density which provides the minimum x 2 value. 
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C. Variation of Nqi/Nqii as a function of Go/™ and $pah 



In Figures [T] through H] we show the variation in Nci/Ncn versus /h 2 as a function of 
Go/n and $pah- In the limit of low Go/n the destruction of C + is dominated by recombina- 
tion on PAH~ and rate equation (1F3 13j) applies. From equation (IB15P we have 



Ncu 



PAH 



n 



Go 



(CI) 



We see in Figures [T] through HI along curves of constant G , the Nci/Nqu ratio rises steeply 
as n increases. 

In the limit of high Go/n the destruction of C + is dominated by charge exchange with 
PAH and rate equation ( 1B16I) applies. From equation ( IB 171) the Nci/Nqh ratio is given by 



Nci 
Ncu 



oc $ 



PAH 



n 



Go 



(C2) 



and the curves start to flatten as n decreases. 



At low values of $pah the destruction of C + is dominated by gas phase recombination 
with electrons, n e . In the limit of low values of Go/n, the electron abundance is determined 
mainly by the photoionization and recombination of Carbon. From equation (1B17|) the 
N C i/N cu ratio is given by 

oc — . (C3) 



ci 



en 



Go 



In the limit of high Go/n, the electron abundance is mainly determined by the soft X-ray 
photoionization of H. In this case n e oc n 1 / 2 and 



N C i n 

oc 



1/2 



Ac ii G 

and we see that the curves flatten considerably at low values of n. 



(C4) 
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Table 2. Reaction Rates 



Reaction 




Rate 
(cm 3 s , 




PAH" + H+ 


-> PAH + H 


8.3 x 10" 


~ 7< £pAH?2~ 


-0.5 a 


PAH + H+ - 


-> PAH+ + H 


3.1 x 10" 


$PAH 




PAH+ + e -» 


PAH 


3.5 x 10- 


- 5 $ PAH T 2 - 


-0.5 


PAH + e -> 


PAH- 


1.3 x 10- 


" 6 $PAH 




PAH- + C+ 


-> PAH + C 


2.4 x 10- 


" ?( ^PAH^2~ 


-0.5 b 


PAH + C+ - 


-> PAH+ + C 


8.8 x 10- 


-9<T) b 
^PAH 




PAH +hu- 


-> PAH+ + e 


4.6 x 10- 


- 9 G exp( 


-1.8A0 c 


PAH" + hv - 


-> PAH + e 


1.2 x 10- 


- 8 G exp( 


-1.8A V ) c 


C + hu^C 


f + e 


2.1 x 10- 


- 10 G exp 


[~2.6A Y ) c 


C+ + e -> C 


+ /iz/ 


1.8 x 10' 


-llT-i-0.83 
1 2 





a T 2 = T/100 K. 



Additional collisonal rates scale as (to) where m is the mass of the collision partner. 



c Photo rates have units s 1 . Go is the F UV field measu red in units of the iHabing (119681 ) 
field (= 1.6 x 10~ 4 erg cm -2 s" 1 sr -1 ). The iDraind (119781 ) field is a factor of 1.7 larger. 
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Fig. 1.— Column density ratio Nci/Ncn versus molecular hydrogen fraction /h 2 = 
2A^h 2 /[^Vhi + 2A/h 2 ]. Observations are shown for cloud column densities in the range 
0.03 < A y < 0.25, JV H2 < 10 17 cm" 2 , and molecular fraction / Ha < 10~ 4 . Clouds with 
Ay < 0.1 are indicated with an "x". Curves show model results for cloud column Ay = 0.1, 
and various values of $ P ah/-R-i7 where R-n = R/10~ 17 . $pah = 2.0, R = 2 x 10 -17 cm 3 
s -1 ($pah/-R = 1; dot), $pah = 1, R — 1 x 10~ 17 cm 3 s -1 ($pah/-R = 1; dash- dot- dot), 
$ PAH = 0.5, i2 = 1 X 10~ 17 cm 3 s -1 ($pah/# = 1/2; long dash), $ PAH = I, R = 2 x 10~ 17 
cm 3 s" 1 ($ PA h/-R = 1/2; dash-dot), $ PAH = 1, R = 3 x 10" 17 cm 3 s" 1 ($ PA H/i? = 1/3; 
short dash)), $ PA h = 0.5, R = 3 x 10~ 17 cm 3 s -1 ($ PA h/-R = V 6 ; solid), and $ PA h = 0.01, 
R = 3 x 10~ 17 cm 3 s -1 ($pah/-R — 1/300; dot; lower curves). The ratio Gq/u varies along 
each model curve with higher Go/n yielding smaller values of /h 2 . Individual models are 
shown with n = 10, 20, 30, and 40 cm -3 (with the exception of the dash-dot and dash-dot-dot 
curves where n = 5. 10. and 20 cm -3 ") and Gr> = 1.7 (O: Blue) and Gr> = 5.1 (A: Red) . 
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Fig. 2. — Column density ratio Nqi/Ncii versus molecular hydrogen fraction /h 2 = 
2Nh 2 /[Nyh + 2Nh 2 ]. Observations are shown for cloud column densities in the range 
0.25 < Ay < 0.75, and N B2 > 10 18 cm -2 . Clouds with A y < 0.5 are indicated with an 
"x". Curves show model results for cloud column Ay = 0.5, R = 3 x 10 -17 cm 3 s _1 , and 
^pah = 1-0 (dash), $pah = 0.5 (solid), and $pah = 0.01 (dot). The ratio Go/n varies along 
each model curve with higher Go/n yielding smaller values of /h 2 . Individual models are 
shown with n = 10, 30, 100, and 300 cm' 3 and G = 1.7 (O; Blue) and G = 5.1 (A; Red). 
For f U2 = 0.02, Go/n ~ 1.7/10 and for / Ha = 0.2, Go/n ~ 1.7/30. 
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Fig. 3. — Column density ratio iV CI /iV C n versus molecular hydrogen fraction /h 2 = 
2Nn 2 /[Nm + 2Nn 2 ]. Observations are shown for cloud column densities in the range 
0.75 < Ay < 1.25 and N Hi > 10 18 cm" 2 . Clouds with Ay < 1.0 are indicated with an 
"x". Curves show model results for cloud column Ay = 1.0, R = 3 x 10~ 17 cm 3 s , and 
$pah = 1-0 (dash), $pah = 0.5 (solid), and $pah = 0.01 (dot). The ratio G /n varies along 
each model curve with higher Gq/ti yielding smaller values of /h 2 - Individual models are 
shown with n = 10, 30, 100, and 300 cm" 3 and G = 1.7 (O; Blue) and G = 5.1 (A; Red) . 



-32 - 




Fig. 4. — Column density ratio iV CI /iV C n versus molecular hydrogen fraction /h 2 = 
2Nji 2 /[Nm + 2Nn 2 ]. Observations are shown for cloud column densities in the range 
1.25 < Ay < 2.13 and N Hi > 10 18 cm" 2 . Clouds with Ay < 1.5 are indicated with an 
"x". Curves show model results for cloud column Ay = 1.5, R = 3 x 1CT 17 cm 3 s , and 
$pah = 1-0 (dash), $pah = 0.5 (solid), and $pah = 0.01 (dot). The ratio G /n varies along 
each model curve with higher Gq/ti yielding smaller values of /h 2 - Individual models are 
shown with n = 10, 30, 100, and 300 cm" 3 and G = 1.7 (O; Blue) and G = 5.1 (A; Red) . 
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Fig. 5. — x 2 pl°t °f ^pah versus H 2 formation rate R. Contour levels are x 2 — 2, 3, 4, 5, 6, 8, 
and 10. Observations are restricted to cloud column densities in the range 0.03 < Ay < 0.25 
and A^h 2 < 10 17 cm -2 . The minimum value of Xmm = 1-6 with 10 sources included. 
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Fig. 6. — x 2 pl°t °f ^pah versus H2 formation rate R. Contour levels are \ 2 = 1.5, 2, 
3, 4, 5, 6, 8, and 10. Observations are restricted to cloud column densities in the range 
0.25 < A v < 0.75 and JV H , > 10 18 cm" 2 . 
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Fig. 7. — x 2 pl°t °f ^pah versus H 2 formation rate R. Contour levels are x 2 — 1) 1-5, 2, 
3, 4, 5, 6, 8, and 10. Observations are restricted to cloud column densities in the range 
0.75 < A v < 1.25 and iV Ho > 10 18 cm~ 2 . 
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Fig. 8. — x 2 pl°t of $pah versus H 2 formation rate R. Contour levels are \ 2 — 1-6, 2, 
3, 4, 5, 6, 8, and 10. Observations are restricted to cloud column densities in the range 
0.25 < Ay < 2.13 and Nn 2 > 10 18 cm -2 . The minimum trough is given by $pah = 
0.22(i?/10- 17 cm 3 s" 1 ) - 0.40 
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Fig. 9. — Mean density plot of $pah versus H2 formation rate R. Contour levels are < n >= 
5, 10, 15, 20, and 30 cm -3 . Observations are restricted to cloud column densities in the 
range 0.03 < A y < 0.25 and iV < 10 17 cm" 2 . 
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Fig. 10. — Mean density plot of $pah versus H2 formation rate R. Contour levels are 
< n >= 15, 20, 30, 40, 50, 60, 70, 80, 90, and 100 cm -3 . Observations are restricted to cloud 
column densities in the range 0.25 < Ay < 2.13 and Nr 2 > 10 18 . The minimum trough from 
Fig. flS]) is given by $ PAH = 0.22(i?/10- 17 cm 3 s" 1 ) - 0.40 
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Fig. 11. — y 2 plot of_$PAH versus H2 formation rate R using C + recombination rate from 



Nahar fc Pradhanl (119971 ). Contour levels are x' 2 = 1.7, 2, 3, 4, 5, 6, 8, and 10. Observations 



are restricted to cloud column densities in the range 0.25 < Ay < 2.13 and N R2 > 10 
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Fig. 12.— Column density ratio Nci/Nqii versus molecular hydrogen fraction /h 2 — 
2iVH 2 /[A^Hi + 2JVh 2 ]. Observations are shown for cloud column densities in the range 
0.03 < A v < 0.25 and N R2 < 10 17 cm" 2 . Clouds with Ay < 0.1 are indicated with an 
"x". Curves show model results for cloud column Ay = 0.1, and $pah = 0.5, R = 3 x 10~ 17 
cm 3 s -1 ($pah/-R =1/6; solid). The ratio G /n varies along each model curve with higher 
G /n yielding smaller values of /h 2 - Individual models are shown with n = 10, 20, 30, and 
40 cm" 3 and G — 5.1 (A) . Arrows show the affects of lines of sight which graze a larger 
(Ay = 1) cloud ('G'), have cosmic rays enhanced by a factor of 10 ('C'), and have half the 
column density in the WNM phase ('W'). 



